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Abstract: 

Usually, Hawking radiation is derived assuming (i) that a future eternal event horizon 
forms, and (ii) that the subsequent exterior geometry is static. However, one may 
be interested in either considering quasi-black holes (objects in an ever-lasting state 
of approach to horizon formation, but never quite forming one), where (i) fails, 
or, following the evolution of a black hole during evaporation, where (ii) fails. We 
shall verify that as long as one has an approximately exponential relation between 
the affine parameters on the null generators of past and future null infinity, then 
subject to a suitable adiabatic condition being satisfied, a Planck-distributed flux of 
Hawking-like radiation will occur. This happens both for the case of an evaporating 
black hole, as well as for the more dramatic case of a collapsing object for which 
no horizon has yet formed (or even will ever form). In this article we shall cast the 
previous statement in a more precise and quantitative form, and subsequently provide 
several explicit calculations to show how the time- dependent Bogoliubov coefficients 
can be calculated. 
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1. Introduction 



Ever since Hawking's original 1974 derivation that black holes emit a steady Planck- 
distributed flux of quanta |2|, there has been a steady and continual stream of 
articles that re-derive the Hawking flux in various different ways |3], £|, ||, [7[ 
|, [TTJ 0, 0, 0, [lj, f§ |T3, 0, 0, [20|, 0, g 0, |2§ — the justification 



being that in doing so one might strip the calculation to its essence and so discover 
what aspects of black hole physics are truly important for the phenomenon, and what 
aspects can safely be put aside. Several important points are by now well established, 
though often not well appreciated: 

• Hawking radiation is more ubiquitous than Bekenstein entropy, and Hawking 
radiation will still occur in situations where the notion of black hole entropy 



has no meaning [10[ [11], [1 



Eternal black holes, and their associated bifurcate Killing horizons, are useful 
mathematical models |3|], and the source for many useful heuristics p6| , but 
physically they do not accurately reflect the formation and evolution of real 
astrophysical black holes. 

The subtle differences between various forms of horizon present in general rel- 
ativity (event, apparent, isolated, dynamical, trapping, etc..) are physically 

and precisely which (if any) of these is 



important @,pL 27, 28 



, B H 

essential for Hawking radiation has direct impact, for instance on the question 
of "information loss" [B3|, |34L B5L BBL B7|. 



In the current article we will address two specific questions (throughout this pa- 
per we will adhere to strictly general relativistic analyses, not considering modified 
dispersion relations): 

(i) How can we derive the Hawking flux emitted by a slowly evolving (as opposed 
to static) black hole-like object? 

(ii) What type of horizon (if any) is required to generate a Hawking-like Planck- 
distributed flux of quanta? 

We shall do so by adapting, modifying, and extending the by now reasonably well 
established result that if (in spherical symmetry) you have an (approximately) ex- 
ponential relation between the affine parameters U and u on the null generators of 
past and future null infinity, J 2 ^ and J^ + , then Hawking radiation will occur. The 
standard form of this argument (as per Hawking's original calculation |], 0) is to 
say that if asymptotically 



U w U„ - A e 



XI) 
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as u — > +00 for some arbitrary positive constants A and Kh, then Hawking radiation 
happens, with fixed non-evolving temperature 

k B T H = ^* (1.2) 

Here U and u can be viewed as different labels that we attach to a null curve con- 
necting J*~ with J^ + . We can equivalently write 

Of course the original Hawking calculation [p], § explicitly assumes the existence of 
an event horizon at Ujj, and that the resulting black hole quickly settles down to a 
static configuration. 

In counterpoint, in 1987 Hajicek || demonstrated that a strict event horizon was 
not necessary, and that a (suitably long-lived) apparent horizon was quite sufficient 
to generate a Hawking flux. (See also [[12], |33], 31].) More recently (within the con- 



text of "analogue spacetimes" [H], [38], ]3j|) the present authors have demonstrated 
that apparent/ trapping horizons can also be dispensed with, or their appearance 
postponed indefinitely into the future E9]. We shall now make these results more 
general, precise, and quantitative pD|| . Inspired in particular by the work of Hu [41]], 



we focus on the existence of an (in our case, approximate) exponential relation be- 
tween the affine parameters on past and future null infinities as the necessary and 
sufficient condition for generating a Hawking flux. 

Here is a summary of our key result: Consider null curves starting from J?~ and 
arriving on There will be some relation between the affine parameters U on 

,JP~ and u on J^ + : 

U = p(u); u = p' 1 (U). (1.4) 

Now pick a particular null curve, labelled by £7* on and on J^ + . We can 
without loss of generality write 



U = U* + C* 







/ ex p 


— I k(u) alu 


1 u* 


J u, 



alu, (1.5) 



for some constant C* and the function k(u) = —p{u)/p{u). Assume (and this is 
where the physics comes in) that k(u) satisfies an "adiabatic condition" 

|«(«*)l < k(m*) 2 . (1.6) 

Then we shall show that this is sufficient to guarantee (under mild technical assump- 
tions) the existence of a Hawking-like Planck-spectrum of outgoing particles reaching 
at u*, now with a time- dependent Hawking temperature 

feflW = ^r^- (1-7) 

Z7T 
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As we move along (that is, as w* increases, possibly even with w* — > +00), this 
relation continues to hold, with the Hawking temperature controlled by as long 

as the adiabatic condition continues to hold. We do not need to assume a horizon - 
of any sort — ever forms. The rest of this article will be devoted to a detailed proof of 
this result. Along the way we shall revisit (and hopefully clarify) the salient features 
that go into calculating the relevant Bogoliubov coefficients — specifically we shall 
very carefully look at the issue of defining and accurately estimating appropriate 
time- dependent Bogoliubov coefficients. 



2. Motivating the adiabatic condition 

Why might we even expect something like the "adiabatic condition" ( |1.6| ) to either be 
true or relevant? We can physically interpret the "adiabaticity condition" as equiv- 
alent to the statement that a photon emitted near the peak of the Planck spectrum, 
with hujao ~ ksTn, that is ~ k, should not see a large fractional change in the 
peak energy of the spectrum over one oscillation of the electromagnetic field. (That 
is, the change in spacetime geometry is adiabatically slow as seen by a photon near 
the peak of the Hawking spectrum.) It is this slow change in the spacetime geometry 
that ultimately permits us to apply a variant of Hawking's original calculation. 

To then verify that this adiabatic condition holds for macroscopic black holes 
let us (for example) think of the standard Hawking calculation for a Schwarzschild 
black hole. The Hawking temperature is |T|. |2|: 

M 2 

k B T H ~hK~-^, (2.1) 

and consequently {assuming self-consistent back-reaction and that the Hawking for- 
mula continues to be true for a slowly evolving almost-Schwarzschild black hole) the 
mass loss rate for a black hole evaporating into vacuum is given by the standard 
result 

A/f 4 

| M | ~ (2.2) 

In particular, as the black hole evolves its temperature changes. This brings up 
and reinforces an important point: Any truly fundamental derivation of the Hawk- 
ing effect should be able to deal with a time-dependent Hawking temperature. If 
your favourite derivation is intrinsically incapable of dealing with time dependent 
situations, then such a derivation is missing fundamental parts of the physics. 

Looking at the surface gravity of the Schwarzschild black hole we can estimate 

K ^Planck f Q \ 

~ ~M^- (2 ' 3) 
So the standard Hawking process for standard Schwarzschild black holes does satisfy 
the "adiabaticity condition" we have enunciated above, at least as long as the black 
hole is heavier than a few Planck masses. 
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Overall this now provides us with a coherent physical picture all the way down 
to the Planck mass, where we see that adiabaticity breaks down, and "quantum 
gravity" (in the sense of "that quantum theory that approximately reduces to general 
relativity in an appropriate limit") takes over. 



3. The exponential approximation 

3.1 Definitions and exact results 

Consider an asymptotically flat spherically symmetric spacetime with a Minkowskian 
structure in the asymptotic past. (The discussion that follows applies equally well to 
any number of spatial dimensions and can easily be generalized to deal with acoustic 
spacetimes in 1+1 dimensions having two asymptotic regions 

|2|.) In the {t^} 



sector of the geometry we define an affine parameter W on J?~ , and use it to label 
the null curves travelling towards the centre of the body. Similarly, u is taken to 
be an affine parameter on J^ + , used to label the null rays travelling away from the 
central body. The independent coordinates {W, u] provide a double-null cover of the 
relevant parts of spacetime (the domain of outer communication). 

As is standard, one can define a canonical functional relationship connecting 
with by using null curves that reflect off the centre at r = 0. This relation can 
be expressed as 

U = p(u), u = p-\U), (3.1) 

where the labels {U, u} are now no longer to be thought of as independent coordinates 
but, since we have explicitly linked them via the function p(-), as different ways of 
labelling the same null curve once it is reflected through the origin. It is to be 
understood that p^ 1 ^) need not be defined on all of J?~ if a true event horizon 
indeed forms; however this function will certainly be well defined on those parts 
of that lie in the domain of outer communication. We shall soon see that the 
function p(-), or equivalently its inverse, is sufficient to encode all the relevant physics 
of Hawking radiation. Specifically, let us choose a reference null curve completely 
traversing the body. It is labelled by on its way out of the body, and by [/* on its 
way in. We want to use "local" information from the vicinity of this reference null 
curve to study Hawking-like radiation that reaches in the vicinity of w*. 

Let us now start the technical computation by defining a quantity k(u) via the 
relation 

K(B) = _;fw = .m (3.2) 

au/au p[u) 

so that k(u) is simply a functional parameterization of the relationship between the 
affine parameters U and u. When this function happens to be almost constant it 
controls the e-folding relationship between u and U, and so provides a notion of 
"surface gravity" in terms of the "peeling" properties of null geodesies. (And so is 
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intimately related to ^effective as defined in |^8, |29| . See also [|3[ [|3, |5| for comments 
on the importance of these "peeling" properties.) Now pick some generic null curve 
labelled by w*, then through integration one can express any U = p(u) as a function 
of its corresponding n(u): 



u = u* + c* 



exp 



k(u) du 



du, 



(3.3) 



for some constant C*. Note that it is impossible, even in principle, to use local physics 
to specify a unique normalization for C*. This is ultimately due to the fact that any 
constant multiple of an affine null parameter is still an affine null parameter. Note 
also that C* is a constant in the sense that it is the same for all null curves u in the 
vicinity of the null curve specified by it*. However C* does depend on the choice of 
the specific null curve one is working around. This formalism continues to make 
perfectly good sense even for k = 0, where it simply implies a linear relation between 
u and U. Let us now write 



k{u) = k(u^) + 8k(u) = + 8k(u)] 
then in particular we have 

r-u 

k(u) du = k*(u — u*) + / 5k(u) du, 
so that (still an exact result) 



(3.4) 



(3.5) 



exp 



k(u) du 



du 



exp [—k*(u - «*)] exp 



8k(u) du 



du. (3.6) 



Here we are interested in situations in which the second exponential on the RHS is 
in some suitable sense (to be more carefully defined below) "close to unity" . 

3.2 Introducing the approximation 

We now begin the approximation procedure: Let us now suppose that 



5k(u) du 



< e 2 < 1, 



(3.7) 



where we shall soon check the conditions under which this happens. Under this 
hypothesis we can re-write the exact result (p.6|) as 



exp 



a* 



k(u) du 



du 



exp [ — k*{u — u*)] {l + 0(e 2 )} du, 



(3.8) 



which we can integrate to yield 



The analysis here is somewhat delicate, because one is integrating a small quantity 
over what could be a very long time. Note that the way we have set things up, this 
approximation will always be valid over some interval — the only real question is 



how long this validity interval will be. (See section 3J3 below.) 

The net result of the discussion up to this point is that, if we accept the condition 
(P-7|), then for u sufficiently close to it* we can effectively replace k(u) by k*, and so 
write 

ru 

U « [/* + C* / exp [— (u - u*)] du 

J u* 

C 

— U* {exp [ — k* (u - u*)] - 1} 

ft* 

[/* + — } - (— e K *"*l exp (-«,«) 
= U* H — A„ exp (—k*u) , (3.10) 

where we have defined 

U* H = u* + — and A * = — e K * n *. (3.11) 

In spite of the similarity with Hawking's approximation, it is vitally important to 
note that Ujj is not the location of the horizon (extrapolated back to ^?~) — it is 
instead the best estimate (based on what you can see locally at w*) of where a horizon 
might be likely to form if the relation between U and u keeps e-folding in the way 
it is at it*. There is no actual implication that a strict event horizon (or indeed any 
sort of horizon) ever forms, only that it "looks like" a horizon might form in the "not 
too distant future". Once we have this approximate relation, 

U = p{u) w U* H - A* exp (-«„u) , (3.12) 

which we shall refer to as the "exponential approximation", then the rest of the 
calculation simply drops out (in the quite usual manner). The only tricky point lies 
in estimating the range of validity of this "exponential approximation" . 

3.3 Range of validity of the approximation 

The exponential approximation condition (|3.7|) , can always be satisfied for small- 
enough integration intervals (u — w*). But how small is small-enough? As a mild 
technical assumption, let us consider only functions k such that we can define a 
constant 

1 l/(n+l) 

D < +oo. (3.13) 



D := sup 



n>0 



1 I ft/* I 



(n + l)\ k^ +1 



Physically this amounts to the assumption that the only two relevant (reciprocal) 
timescales in the problem are and Dk*. (Any other scale is assumed to be smaller, 
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which is simply another way of saying that k(u) is slowly varying over the region 
of interest.) This condition is rather mild, covering even functions with poles at 
specific values of it. So k(u) is even allowed to exhibit "sudden singularities", with 



this terminology being borrowed from cosmology [26, J7, 48 
Under this hypothesis we have 



pu 




I 5k(u) du 









+oo 1 

^(n+l)! * 

71=1 v ' 



+oo 



\u — u, 



n=l v ' 
+oo 

< D n+1 K™ +1 \u — 



in+1 



(3.14) 



71=1 



Let us temporarily set x = Dn*\u — u*\, then 



Sk(u) du 



Now, as long as 



<5>' 

n=2 



X 



1 — X 



1 — X 



X 



1 — X 



1 — X 



< e l 



(3.15) 



(3.16) 



we are sure that condition fl3.7| ) is satisfied. Taking into account that e <C 1, the 
previous condition is certainly guaranteed to hold as long as x 2 < e 2 /2. (This is 
not the optimal condition, but it is simple and quite good enough for our purposes.) 
Then 

2D 2 kI (u - u*) 2 < e 2 < 1. (3.17) 

Thus the range of validity of the exponential approximation condition ( |3.7|) is cer- 
tainly at least as large as 

el , 

\u-uA< — < ^= • (3.18) 

In the most simple situations the first term in the definition (|3.13|) dominates, thereby 
yielding 

(3.19) 



2D 2 



1 ' 



Hi 



For instance, for an evaporating Schwarzschild black hole one can estimate 

M 4n+2 



«(»)| ~ l X 4 x 7 x ■ ■ ■ x (3n - 2) 



Planck 

M 3n+1 ' 



(3.20) 



and so 
1 



(n + 1)! «J 



i <»i 

| | 

l 7^ 



l/(n+l) 



n 



1 M 2 ™ 

J m Planck 

1 M 2 ™ 



l/(n+l) 



M 2 

1V1 Planck 
M 2 



M 



1 2/(n+l) 



Ml 



Planck 



(3.21) 
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It is the trailing factor of (M/Mpi anc k) 2 ^ n+1 ^ that guarantees that for M 3> Mpi anc k 
the lowest derivative dominates in the definition of D. In fact careful estimates 
putting in all relevant numerical pre-factors show that for a Schwarzschild black hole 
the first term (the k term) dominates down to some 5 Planck times before the final 
evaporation event. 

In typical situations like the above, where the k term dominates in the definition 
of D, the exponential approximation holds for an interval at least as large as 



\u 



«*| < t\k,\' 1/2 < \k*\~ 1/2 . (3.22) 



Hereafter we will assume that this is the case. (Other more complicated situations 
can be dealt with in a tedious but straightforward manner.) 

What is the physics behind this condition? 

• If one looks at the spacetime for a limited interval Au = \u — u*\ then the 
lowest possible frequency one can hope to resolve is fi min ~ 1/Au. 

• The condition \k*\ (u — u*) 2 < e 2 <C 1 then collapses to the statement that over 
the time interval of interest A/t < e 2 f2 min <C fi min . 

• That is, the "exponential approximation" is guaranteed to hold as long as the 
change in surface gravity is less than the minimum detectable frequency. 

Over what range of situations does this condition hold for macroscopic black holes? 
For Schwarzschild black holes our previous estimate of \M\ implies 

/If 6 

k _ ^Planck (3>23) 

so the condition for guaranteed validity of the exponential approximation reduces to 

/ M \ 2 f M \ 2 

\u - u„| < t Planck x — ~ 10 15 x (age of universe) x — . (3.24) 

V J- VI Planck / \ M & J 

So for stellar mass black holes there is absolutely no difficulty in satisfying this 
inequality. 

Note that we can rewrite the condition 

|k*| iu-u,) 2 < e 2 < 1, (3.25) 



as 



[k^u-u^Y < e 2 < 1. (3.26) 



Since we do want to be able to have /t* \u — uj\ ^> 1, (that is, as we have just seen, 
macroscopic mass black holes should be analyzable via quasi-static techniques for 
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many light crossing times), the condition for validity of the exponential approxima- 
tion in turn implies 

< e 2 ; whence « 1. (3.27) 

Note now \k*/ nl\ is "very very much less than 1". In this sense validity of the 

exponential approximation implies the adiabaticity condition. 

It is now useful to define an interval S + in terms of the affine parameter on J ?r+ , 



S 4 



( 



-1/2 



6 



-1/2 



)■ 



(3.28) 



and the corresponding interval <S_ defined in terms of the affine parameter on , 
S_ = p(S + ) = (u* H - — exp [e«*J/«*| _1 / 2 ] , Uff - — exp [— e/c*|/c*| _1 / 2 ] J . 

\ ft* ft* J 

(3.29) 

The "exponential approximation" is then valid for u e S + , which is equivalent to 
U £ S-. Note that for a/|k*|/k* sufficiently small we can get "arbitrarily close" to 
the "might be horizon" at Ujj while still remaining within the realm of validity of 
the exponential approximation. 

3.4 Nonlocal normalization 

To determine the constant C*, and so also ultimately determine the constant A* ap- 
pearing in the exponential approximation, requires an overall choice of normalization 
- we need to relate the overall scale of the affine coordinate on to that for J? + . 
A natural choice is to demand 



p(u — > — oo) u. 



(3.30) 



This has the effect of making sure that in the infinite past (i.e., before any collapse 
or other dynamics was initiated) and are connected in the simple sensible 

way 

dU 

— > 1 as u — > — oo. (3.31) 

du 

But from 



dU 

— — = G* exp 
du 



we then see 



exp 



n(u) du 



n(u) du 



(3.32) 



(3.33) 



So we see that C* depends on the entire past history of J? + (the history of the 
collapse in Hawking's language). Furthermore we now estimate that a horizon might 
form at 

i r f u * 1 

Uh{u*) = U* H = U* H exp — / k(u) du , (3.34) 
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and that 



1 

A* = — exp 

ft* 



k(u) du + ft* Hi, 



— oo 



(3.35) 



So finally we have explicit formulae for the parameters appearing in the exponential 
approximation. 

We now also have a clear physical interpretation for C*: Since 



dU 
du 



C*, (3.36) 



tit 



we see that a small time interval AU around U* on J" is related to a small time 
interval Am around m* on ^ + by 

AU = C* Au. (3.37) 

This implies that C* is the Doppler shift factor (redshift factor) relating the frequency 
of a photon emitted at £7* on to its frequency when it reaches w* on J^ + . 

3.5 Logarithmic approximation 



One can also invert the exponential relation (|3.12 ) to give the perhaps more common 
"logarithmic approximation" QTj, 

U ^J_ i n ( U -JLZ^- 1; f/G5_. (3.38) 



For completeness, we note that in terms of the formalism we have developed above, 
the definition 

, , d 2 U/du 2 

"<"> = -5/7*r • (M9) 

can, after some manipulations, be seen to be equivalent to 

. * d 2 u/dll 2 , , . 

K(c,) = ' (3 ' 40) 

which can be integrated as 

= ^ + / f fr - ( 3 - 41 ) 

with the same C* as occurred previously. This can be used as a basis for an alternative 
way of deriving the "logarithmic approximation" above. 
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4. Peeling properties versus inaffinities 



To obtain a clear physical interpretation for the function k(u), we will now connect 
the behaviour of the affine null parameters u and U back to the geometry as encoded 
in the spacetime metric. The most efficient way of doing so is by adopting double null 
coordinates {W, u] on the transverse slice (the {t,r} sector). We use W as the null 
coordinate rather than U, to emphasise that W and u are to be viewed as independent 
coordinates on the spacetime, whereas U = p(u) and u are linked by being different 
labels for the same null curve (reflected through the origin). We wish specifically to 
clarify the relationship between k(u) and the "peeling properties" of outgoing null 
geodesies, and contrast this with the inaffinity properties of null geodesies on the 
horizon. See the appendix for additional background. (See also |42| and [[43| for 



comments on the importance of these "peeling" properties in "analogue spacetimes, 



and pi , 4"5| for comments in a general relativistic context. 



4.1 Metric asymptotics 

For any spherically symmetric spacetime, we can without loss of generality adopt 
such double-null coordinates and set 

ds 2 = -F(W, u) dW du + r(W, uf dn 2 d _ v (4.1) 



See, for instance, Bergmann and Roman f2"? |. Here the two functions F(W,u) and 



r(W, u) completely specify the spacetime geometry. Since we have chosen W to be 
affine on ,J?~ and u to be affine on we can impose 

lim F(W,u) = l; lim F{W,u) = 1. (4.2) 

W— >+oo u — y — oo 

This makes the metric particularly simple on J? . 
4.2 Inaffinity estimation in the bulk 

Recall that the generators of are null geodesies, affinely parameterized by W. 
Using (W, u, 6*i, ... , 6d-i) coordinates, throughout the spacetime one can define the 
null vector field 

fc a = (1,0,0, •■■,()), (4.3) 

pointing in the direction of increasing W . For that null geodesic that skims along 
this is an affine parameter, for other null curves (labelled by a constant value of 
u) which traverse the interior of the spacetime this is no longer an affine parameter. 
The 4-acceleration 

k a V a k b = T b ww = ^(l,0,---,0), (4.4) 
now enables us to identify a new and logically distinct quantity 

«buik(W,u) = % ; . (4.5) 
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This is now a "bulk" quantity denned everywhere in the spacetime, (not just on 
horizons or asymptotic infinity). This Kbuik is a "bulk inaffinity estimator" which 
measures the extent to which the coordinate W fails to be an affine parameter along 
null geodesies of increasing W. It is this "bulk inaffinity estimator" that is closely 
related to the textbook notion of surface gravity, while it is the "peeling" notion of 
surface gravity, k(u), that we have seen is related to the Hawking flux. 

4.3 The surface gravity of T-L + 

If (and this is one of Hawking's key assumptions in his derivation of the existence of 
a Planck flux reaching J^ + ) a true future-eternal event horizon forms, then there will 
be a region of J?~ (namely W > Uh) which has no natural "lift" to <J? + . Instead, 
this region of "lifts" to "H + , the future-eternal event horizon. On "H + we have 
an (in principle evolving H^-dependent) inaffinity notion of surface gravity: 

KWEnityW = lim K hu i k (W,u). (4.6) 

u— »+oo 

Note particularly that by construction the textbook notion of surface gravity makes 
sense only on the future event horizon % + . In contrast, note that our definition k{u) 
naturally resides on not on T-L + . 

It is only if we make the further assumption that, after the future-eternal event 
horizon forms, the black hole settles down to an asymptotically static state, that we 
have a relation 

k h = lim ftinaffinityW = lim k(u), (4.7) 

W— ¥+oo u— >+oo 

where kh is now the quantity appearing in Hawking's original derivation. We em- 
phasise that without such an asymptotic assumption the two notions (inaffinity vs 
peeling) are distinct, and in general different. Even in this particular case they at 
best only coincide at i + , future timelike infinity, and are unrelated at other locations. 
It is only for static (or stationary) spacetimes that the two notions exactly coincide 
for all times. (See appendix \K72^ ). 

If in contrast one has a future non-eternal event horizon, (an event horizon that 
forms and subsequently evaporates as described by the standard conformal diagram 
for an evaporating black hole, see for instance |13|), then as seen from ^ + the horizon 



both forms and evaporates "simultaneously" at some uh < +00, and one has 

U H = p{u H ), U E =p(u + H ) } (4.8) 

where Ue now labels the first incoming null ray that reaches the centre after the 
evaporation process is complete. Note that at the end-point % of the evaporation 
p{u) has a discontinuity, and n{u) will diverge. This strongly suggests that in this 
situation some form of "thunderbolt" will be emitted ffM FOI . 
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4.4 Summary 



The physics here is intriguing |4(| - to get a Planck flux not only do we not ever 
need the future horizon T-L + to form, but the Hawking temperature is not logically or 
physically connected to the surface gravity of the horizon ~H + . Instead the Hawking 
temperature is primarily related to k(u), that is, the "peeling off" properties of the 
null geodesies that actually do reach J? + , and hence to the peeling notion of surface 
gravity. It is only when a collapse settles down to a static black hole spacetime (the 
end point of classical collapse in general relativity) that these quantities happen to 
be equal. 

5. Bogoliubov coefficients: Basic framework 

The Bogoliubov coefficient calculation — which proceeds from the exponential ap- 
proximation to the Hawking flux — is in principle completely standard [|TJ 0, |5j]] . 
Modulo normalizations and technical issues, such a calculation will automatically give 
us a steady Planck flux of particles at as long as the exponential approximation 
holds for a sufficiently large interval. The purpose of presenting the calculation in 
some detail is to provide a very general and clean derivation, exhibit some novel re- 
sults using the stationary phase technique for evaluating the relevant integrals, and 
to pull together a few tricks otherwise scattered throughout the literature. 

To place this discussion in a textbook perspective: Note that the ultimate goal 
is to reproduce some close analogue of the discussion in Birrell & Davies [fyj ; either 



equations (4.59), (4.60), and (4.61) on page 108 of the "moving mirror" discussion 
of pages 102-109, or the closely related equation (8.32) in the collapsing black hole 
discussion on pages 250-260, but now with slowly varying time-dependent Bogoliubov 
coefficients. 

5.1 Bogoliubov coefficients in terms of the Klein— Gordon inner product 

It is very easy to write down the asymptotic modes for a massless scalar field. On 
the properly normalized modes are (see for instance Birrell & Davies equation 
(8.29) and (8.30) on page 258; note we are now in 3+1 dimensions) 

exp(— iujU). (5.1) 



These modes extend to the full spacetime as 4> u (ou; t, x). Similarly, on the prop- 
erly normalized asymptotic modes are 

exp(— icon), (5.2) 



which extend to the full spacetime as <fr u (uj]t,x). The exact Bogolubov coefficients 
are now easily written down (in terms of the Lorentz-invariant Klein-Gordon inner 
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product, see Birrell & Davies equation (2.9) on page 11, and equation (3.36) on page 
46) as 

a(u, J) = (<p u (uj; t, x), u (lu'; t, x)) (5.3) 
= -i I d 3 x (0 u (w;t,x) d t (p u (uj';t,xy - <p u (u';t, x)* d t u (cu;t, x)} , 



f3(u,u') = -(4> u (uj;t,x),4> u (cu';t,x)*) (5.4) 
= -i I d 3 x{(t> u (u';t,x) d t <p u (cu;t,x) - (f> u (u;t,x) d t <f> u {u';t,x)} , 



where the integral runs over any arbitrary spacelike hypersurface that terminates at 
spacelike infinity (that is, at Note that we associate u' with the frequency of the 
mode that is simple on J^~, and oj with the mode that is simple on J^ + . 

Assuming for the moment that no event horizon forms, we now choose the space- 
like hypersurface to skim arbitrarily close to J^ + . (This is a slightly unusual proposal, 
because calculations in the literature, where one commonly assumes that an event 
horizon forms, are normally performed on J>~ |], ^|, |5"lj . However, if no horizon 
forms, the ultimate results are equivalent, because the Bogoliubov coefficient, being 
a Lorentz invariant "inner product" , will not depend on which hypersurface one picks. 
Thus working on is a useful consistency check on the formalism. In particular, 
this way of organizing things emphasises the fact that actual horizon formation is 
not central to the Hawking process.) Focussing on the (3 coefficient we can evaluate: 

^'^=0 — 7rf=i (5 - 5) 

Z7T V 4U It)' 



du < exp[— ioju] d u exp[—iu'p(u)] — exp[— iu'p(u)] d u exp[— iuu] 
Integrate the first term by parts. There will be a surface term proportional to 



exp[—ioju]exp[—iu)'p(u)]\_ 00 , (5.6) 

which can safely be discarded. The remaining term, suppressing the angular Kro- 
necker deltas, gives 



1 fu 



+oo 



(3(u),uj') = — \ — / du {exp[— iuu] exp[— iu'p(u)] } . (5.7) 
2vr V w ./ x 

Compare with Birrell & Davies equation (4.59) on page 108. We emphasize that this 
formula, though it looks extremely simple, is in fact one of the key results needed in 
the derivation of the Hawking flux. We can also perform a change of variables in the 
above, and after an integration by parts obtain the completely equivalent formula 



i fu7 f +oc 

p(u,u/) = — \ — / dU {&v[-%up-\U)\ex$[-iu'U\ } . 
2-k V u y_ 00 



(5i 
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Of course this is also exactly what you would get by evaluating the Bogoliubov 
coefficient by using a spacelike hypersurface that skims just above J?~ . 

If in contrast one assumes that an event horizon does form, then it is certainly 
more efficient to work on . (Since is always a complete Cauchy hypersurface, 
whereas if a horizon forms one needs to work with U "H + to get a complete 
Cauchy hypersurface.) In any case there is universal agreement on the applicability 
of equation (518). See for instance P], % pT]]. 

5.2 Approximation scheme 

The novelty in the current analysis is this: Expand the integrand around U = 
corresponding to m = where we have already seen that we can justify the "expo- 
nential approximation" 

U ^U* H -A, exp(-/c*«). (5.9) 
In view of this, we have the approximate result 

exp[— icu'U] = exp[— iu'p(u)} ~ exp[— iu'U^j] exp {-Mu/A* exp[— k*u]} , (5.10) 

which we know to be valid in the range 

\u-u*\ < e\k*\- 1/2 . (5.11) 

Then expanding the integrand (not the integral) 



/•U*+€|k*| -1 '' 2 

x / du exp[— iuu) exp[+iu'A^ exp{— k*u}]. (5.12) 

Jut.-elkJ- 1 /' 2 



The most significant pieces we are currently neglecting in this approximation are the 
contribution from the "tail" integration regions 

(-00, - e\kJ~ 1/2 ) and (it* + e\k*\~ 1/2 , +oo) . (5.13) 

That is, we are neglecting the contribution from the region M — S + , since in the 
region 

\u - «*| < e\kJ~ 1/2 <^> u e S + (5.14) 

we are accurately estimating the value of the integrand. But in the tail regions the 
phase is rapidly oscillating, so the tail contribution to the Bogoliubov coefficient is 
small. 

Hawking's key approximation |], [2| (suitably rephrased, reinterpreted, and now 
extended to slowly evolving systems) is this: In view of the condition underlying 
the exponential approximation, one is justified in extending this approximation for 
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the integrand to all it, but provided that one then interprets the resulting Bogoliubov 
coefficient as the "effective" /3(it*) for the production of particles (wave packets) that 
arrive at J? + at u ~ u*. Within the context of this approximation we now have 

1 [uj f +oc 
(3(u*;u,u)') rs — \ — exp[—iu}'Uff] / du exp[— iuju + iu'A* exp{— /«*«}]. 
2vr V W > /_ 00 

(5.15) 

One is effectively assuming that the true integrand and approximate integrand are 
both rapidly oscillating in the "tail" integration regions in ( |5.13 ), and so both average 
to zero on these intervals. 

We emphasize that this approximation (though without the extension to time- 
dependent Bogoliubov coefficients) is actually what Hawking implicitly does in his 
original 1974 calculation, with slight modifications since he is working on J?~, and 
for the particular case U^ = Uh, corresponding to formation of a strict event horizon 
where w* — > +oo. In Hawking's Nature article JIJ, the approximation is hidden in 
the third column where he first writes down the approximate asymptotic form of 
the mode and then performs an exact Fourier transform for this asymptotic approxi- 
mation. Hawking makes essentially the same approximation in his Communications 
in Mathematical Physics article ||, where (2.18) is his asymptotic estimate for the 
mode on J?~, and by (2.19) he has performed an exact Fourier transform of this 
approximate asymptotic mode. 

One can also justify the approximation using a wave-packet argument. We know 
the exponential approximation with fixed it* is valid for u G S + . Now consider some 
wave-packet that arrives at J? + with compact support in S + . Physically, such a wave 
packet cannot tell the difference between the (unkown) exact relation u = p~ l {U), the 
"exponential approximation" defined above, and Hawking's (at first glance seemingly 
absurd) extrapolation of the exponential approximation to "all time". That is, for 
any wave-packet with compact support on S + , where by definition the quantity 
k(u) ~ k(u*) is approximately constant, we can use the standard Hawking calculation 
to derive the Bogoliubov coefficients relevant to wave-packets confined to this time 
slice — that is, (3(u*;u,u'). We see that up to an irrelevant overall phase: 

1 fuj f +oc 

/3(it#; tu, u/) ~ — \ — / du exp[— iuju + iu'A* exp{— u}]. (5.16) 
2vr V w y.oo 

At this stage you can see that the standard Hawking result is really unavoidable. To 
complete the calculation we can now either approximate this integral using stationary 
phase techniques, or go for an "exact" evaluation using Gamma functions. 

6. Stationary phase evaluation of the Bogoliubov coefficients 

Stationary phase approximations have been somewhat neglected in this part of the 
technical literature, but appear to be at least as reliable as the usual "exact" results 
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obtained via Gamma function identities. One actually has to do two slightly different 
but very closely related calculations — a Boltzmann calculation to get the overall 
normalization and a calculation of the ratio \/3/a\ to obtain a Planck spectrum. 

6.1 Boltzmann spectrum using stationary phase 

The stationary phase approximation is 



/+oo 
exp[— i(f)(u)]du 
-oo 



2ni 



exp[-i(f)(u s )], 



(6.1) 



where u s is defined by <f>'(u s ) = 0. This represents the leading term in an asymptotic 
approximation based on a Gaussian expansion of the integrand around its stationary 
point. In the current situation the relevant phase is 



4>{u) = uju — u'A* exp(— k*u), 
and the stationary phase condition is given by 

4>'(u s ) = uj + uj'A*k* exp(—K*u s ) = 0. 

That is 



exp(—K*u s ) 



oj 



UJ 1 A* 



whence 



UJ 



"> — — 1,1 % . ,/ j 

\ OJ si^K,* 



Note that u s , the location of the "stationary" phase, is complex. Now 

4>"(u s ) = —u)'A m Klexp(—K m u 8 ) = UJ. 
The stationary phase approximation is now 

1 



(6.2) 

(6.3) 
(6.4) 
(6.5) 

(6.6) 



f3(u*; uj, uu') 



whence 



/?(«*; uj,uj') 



fl7 


2ixi 




V K^UJ 




1 



exp[— iojUg + iuj'A* exp(— k*u s )], (6.7) 



cxp 



UJ 71 
—lUJUs — I h %- 



(6.8) 



ft* 4_ 

The only really interesting bit comes from the imaginary part of u s which yields, up 
to a physically irrelevant overall phase: 



1 



This gives 



/\|2 



(u*;uj,uj)\ 



^J2lXK*UJ 

1 



= exp 



2ttk*u 



7 ex P 



UJ1X 

2tiuj 



(6.9) 



(6.10) 
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which is a Boltzmann factor. Compare the normalization with Birrell & Davies 
equation (4.61) on page 108. Note that this elementary application of the stationary 
phase approximation does not (yet) yield a Planck spectrum, though it does yield a 
Boltzmann spectrum, and one quickly gets to the Hawking temperature this way: 

k B T H (u*) = ^ = (6.H) 

Note that the size of the sub-leading terms in the asymptotic expansion underlying 
the stationary phase approximation will be controlled by the width of the Gaussian 
that is its leading term, that is by <p"(u s ) = uk*. This implies that the stationary 
phase approximation (as we have used it so far) will be valid up to fractional cor- 
rections of order 0(k*/uj). This is why the current version of the stationary phase 
approximation is intrinsically a high-frequency approximation, incapable of seeing 
the sub-leading exponentials critical to probing the Planck nature of the spectrum. 

6.2 Planck spectrum using stationary phase 

To get a Planck (rather than Boltzmann) spectrum use the exact normalization 
condition 

du>' [a(wi, oo')a*{u 2 , J) - ou')f3*{u 2 , u')} = Sfa - u> 2 ), (6.12) 

and consider the ratio 

0(u m ;u,u') f*™ duexpl-iuu + iu'A* exp(-/t*w)] 



a*(u*; u, u') du exp[— icon — iu'A* exp(— k*u)] 



(6.13) 



Now use stationary phase on numerator and denominator separately. We shall see 

that this ratio is much better behaved than each individual integral. The relevant 
phases and derivatives are: 

= urn =)= u'A* exp(— k*u), (6.14) 

0' = uj ± u'A*k* exp(— k*u), (6.15) 

4>" — ^uj'A^kI exp(— k*u). (6.16) 

Furthermore for n > 2 

0(») = T (-l)VA»<exp(-K*«). (6.17) 
The two stationary points differ only in their imaginary part 

u . = -lh,( « )-]S^, (618) 
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which we can rewrite as 



u. 



it* In 



UJ 



u'C* 



+ 



ITi 



For both integrals 
and for n > 2 
while for the phases 



4>"(u s 



Us 



K*UJ, 



4>(u s ) = uu s + lo/k 

Since the two stationary points differ only in their imaginary part, we see that 



(6.19) 

(6.20) 
(6.21) 
(6.22) 



exp(— tvoj/k*), 



(6.23) 



and we are essentially done. Note that the sub-leading terms in the asymptotic 
expansion for j3 are the same as the sub-leading terms in the asymptotic expansion 
for a, so they quietly cancel, and the stationary phase approximation for the ratio 
\/3/a\ is much better behaved than the stationary phase approximation for each term 
individually. In fact, using the normalization found in the previous (Boltzmann) 
version of the stationary phase calculation and showing explicitly their phase factor 
we can write 



/?(«*; uj, U)') 



a [u*;uj,uj 



2ttk*uj' 
1 

exp 

2/KKJjj' 



ft* 


In 




V w ) _ 


. U) 
ft* 


In 







Ik 



(6.24) 
(6.25) 



where a w ,(3 w are real functions of uj such that = exp(— ttu/k*). In this way, 

using C |S. 121 ) it is easy to check that we now get 



A|2 



\f3(u*;uj,uj')\ 



2ttk^lo' exp(27ru;/fi;*) 



(6.26) 



which is indeed a Planck spectrum. Compare, for instance, with Birrell & Davies 
equation (4.61), page 108. 

6.3 Mathematical range of validity for stationary phase 

A tricky point of the stationary phase approximation is this: Is the location u s of 
the stationary point in the phase within the range of validity of the exponential 
approximation? To check this consider 



|Re(w s — u* 



1 



In 



UJ 



Ul' A^K* 



It, 



1 



1* 



111 



Ul 

uTcl 



(6.27) 
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Thus mathematically the stationary phase approximation will be valid as long as 



that is 



exp [— e K 



In 



-1/2 



CO 



u'C* 



-1/2 



]< 



IjJ 



< exp [+e 1 ' 2 ] . 



(6.28) 



(6.29) 



Remember that u/ is the unobserved frequency on J"~ while u; is the physically 
relevant frequency on . Now consider the quantity u'. What is its physical 
interpretation? Consider the null curve labeled by w*, and a small interval Au 
centered on -u*. Then on this interval corresponds, (recall ( |3.37| )), to AU = 
Au, that is 



AH' * 



AU 



Therefore 



(6.30) 



(6.31) 



where u' + is the "Doppler shifted" value of u' if we transport it to y + along the null 
curve labelled by u*. (That is, u' + is the frequency on corresponding to u' on J?~ 
at "time" u*.) Then the condition for validity of the stationary phase approximation 
is 



exp [— 



-1/2 



]< 



^1 



< exp [+6K* I I 1/2 ] 



(6.32) 



This can always be satisfied for an extremely broad range of frequencies. Indeed, 
for a Schwarzschild black hole we have already seen k ~ 1/M and k ~ M| lanck /M 4 , 
therefore 



exp 



eM 



Mi 



Planck 



< 



< 



exp 



eM 



M 



Planck 



(6.33) 



Thus mathematically there is no risk of the validity of the stationary phase approxi- 
mation causing additional problems — the validity of the exponential approximation, 
the adiabatic condition, and the validity of the stationary phase approximation are 
all secure for macroscopic black holes. If we wish to more physically relate u' to u, 
(currently they are completely independent quantities), that can be done once we 
perform a wave-packet analysis. See section l8|. 



7. Gamma function evaluation of the Bogoliubov coefficients 

As an alternative way of evaluating the integral (and hence yet another consistency 
check) we now use Gamma function techniques. Note that the logic to be presented in 
this section is at first sight slightly odd, but completely standard [|], || . We are taking 
an approximate value for the Bogoliubov coefficient and then exactly evaluating the 
resulting integral. The ultimate justification for this comes from looking at a wave- 
packet supported on S + and noting that for this wave-packet one might as well 
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extend the integral to "all time". Discarding an irrelevant overall phase, start with 
the integral 

1 fuj f +0 ° 

/3(u*;uj,uj') ~ — \ — / duexp[— iuu + iu'A* exp(— k*u)]. (7.1) 

Now perform a change of variables 

z = exp(— (7.2) 

so that 

i ATT i 

P(u*;u,uj')n — J-— dzz^l**- 1 exp(ia/A.z). (7.3) 

27T V W Jo 

Now take 

z = -iu/Az, (7.4) 
so after an appropriate change of contour 

i /nr i /" +o ° 

/3(m,;w,w') « — J— — (-iu/A*)-™'*' / dzz™/**- 1 exp(-z), (7.5) 

27T V ^ i 

whence 

w, u') « i/^ — (-iw'A.r^Ttiw/K,). (7.6) 

27T V ^ K * 

So up to an irrelevant phase 

w, w') ~ 7^- — Hr w/K * r(iw/«0 (7.7) 

= — \ —, — exp[— 7tu;/(2k*)] YUu/k^). (7.8) 
27r y w K, 

Using the properties of Gamma functions 

|rfcr)| 2 = .1 (7.9) 



we have 



, 2 1 oj 1 exp(— ttuj/k*) 



so that 



\P(u.;u,u')\ 2 « -— j - - tt " v , ' */ (7.10) 
(27r) 2 uj' k% u/k* sinh(7ra;//c*) 

|/3(u.;u;,u/)| 2 «— * ; r - (7.11) 
2nn*u' exp(27rw//t 4t ) — 1 

Compare with Birrell & Davies equation (4.61), page 108. This is a Planck spectrum 
with 

feW = ^ = ^. (7.12) 
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And we are done. Note however, that this calculation is actually no more "exact" 
than the stationary phase approximation; one has just hidden the approximations 
elsewhere — in particular, an approximate form for the mode function near u* has 
been promoted to all of ^ + . That is: Our approximate mode function should be a 
good approximation to the real mode function for 

\u — u*\ < e|/t*| -1 / 2 . (7-13) 

We then assume this region S + dominates the integral, so that formally we can 
integrate the approximate mode function all the way from — oo to +00. Of course, 
the condition that this region S + dominates the integral is equivalent to requiring 
the validity of the stationary phase approximation, leading again to the constraint: 

to 



exp [— e/c«,|«*| 1 ^ 2 ] < — < exp [+e/c*|«*| 1//2 ] . (7-14) 



8. Physical particle detection: Wave-packets 

When we say that we are calculating the number of particles of frequency u arriving 
to the asymptotic region at time w*, what we really mean by particle is a normalized 
wave packet with a frequency content concentrated around the frequency u while si- 
multaneously being temporally localized around u*. Let us distinguish the frequency 
at which the wave packet is centered, denoting it u, from the dummy label u of 
generic plane-wave-like modes. Being mathematically more precise, what one has to 



calculate is actually a wave-packet representation of the Bogoliubov coefficients (fy] : 



"+00 

AvprO*; w,u;') = / dCoP^.^Co) f}(u*\ Co, u'). (8.1) 



00 



Here P u ,-,uj(w) defines the frequency distribution of the wave packets x), 

/+00 
duP u ^(u)<j) u (u-,t,x). (8.2) 
-00 

These wave packets are taken to be normalized so that 

X) 

\P u ^{uo)\ 2 duo = 1. (8.3) 



In the vicinity of u the function P u *,u{uj) can always be approximated by a slowly 
varying absolute value multiplied by a rapidly varying phase of the form e tU)U *. For 
instance, take a simple wave packet representation 

Pu,A^) = 7— [0(w - oj + Acu/2) -Q(p-u- Acu/2)] e iQu * , (8.4) 
with Ao; the width of the wave packet. 
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Now, in all the previous calculations we were neglecting an irrelevant phase factor 
in u), a/). It is easy to check that this phase factor always contains a term of 
the form 

exp (^i^\n[uj' A^k^/uj]^ . (8.5) 

Taking this into account, it is not difficult to see that for u ^> Au an approximate 
evaluation (the slowly varying factors that depend on u are approximated by their 
value at u) of /3wpr(w*; u>, uj') yields 

sin <^ Au (w* - i \n[u'A*K*/uj] J > 



Aw ^u* — iln[u;M.*K*/u>] 



1 v ' ' 71 (-[Aw/k»] ln[w'C»/a;]) v ; 

Then the number of particles (wave-packets) of frequency u arriving at infinity at 
time is 



/ 



du'\[3 WPR (u*] u,u')\ 2 - 



exp(27rw/K*) — 1 

1 f + °° duj' sin 2 (-[Aw/^] ln[w'C*/w]) 



x 



2itk* J_ 00 u' (-[Aw/k 4 ] ln[o;'a/w]) 

' (8-7) 



exp(27ra;/fi;*) — 1 

which is finite and has precisely the anticipated Planck shape. 

Furthermore, the function sin 2 (x)/a; 2 is tightly peaked around x — 0, with a 
"width" characterized by the location of the first zero at |x| = tt. So the range of 
"important" values of u>' contributing to this integral is characterized by 

[Au/k,] \Hu'CJcu)\<n. (8.8) 

That is 

OJ 

For physically reasonable wave-packets, with physically reasonable widths Au, this 
interval is well inside the interval of mathematical validity of the stationary phase 
calculation as estimated in ( |6.28| ), Q6.29Q , so in (|8.6|) it is perfectly justified to use 



a (3(u*;u,u') calculated through this approximation technique. Note that for wave- 
packets the frequency u' on J"~ is related to the central frequency u on ^ + by 
the redshift factor C* and an envelope depending on the frequency width Aw of the 
wave-packet. 
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9. Physical necessity of the adiabatic condition 



We have now, subject to the approximation underlying the exponential approxi- 
mation, which we emphasize is always valid on a small enough region, derived the 
existence of a (slowly varying) flux of particles arriving at with a characteris- 
tic temperature ks Th{u) = H k,(u)/2tt. The derivation was ultimately based on 
wave-packets of support confined to the region S + over which the exponential ap- 
proximation is valid. But a wave-packet supported on S + is bandwidth limited to 
only contain a limited range of frequencies 



Ac > y^. (9.i) 

If we wish the wave-packet to contain frequencies capable of probing the Planck peak 
of the Hawking spectrum (instead of merely probing the high-frequency Boltzmann 
tail) then we must have Au < «*. Then: 



V 



^<e 2 ; (9.2) 



k-2 ~ ' K-2 



So here we have another viewpoint on the necessity of the adiabatic condition: While 
by direct calculation the adiabatic condition is seen to be automatically valid for 
macroscopic (Schwarzschild) black holes, the reason the condition is useful is that it 
permits the calculation to probe the Planck peak of the spectrum. 



10. Discussion 



Note what we think is the importance of these results — the formalism we have 
developed is an extremely general way of getting Hawking-like fluxes from a rather 
broad class of spacetimes — of course it agrees with previous results and also puts 
things into a more general framework. The key result is this: 

Whenever and J? + are connected by a function U = p(u) of the form 



:io.i} 







exp 


— 1 k(u) du 


lu* 


J u* 



du, 

for some constant C*, and where k(u) satisfies an "adiabatic condition' : 

dn(u) 



du 



< k(u*) , 



(10.2) 



then the argument sketched above proves the existence of a Hawking-like Planck- 
distributed spectrum of outgoing particles reaching at with a (time-dependent) 
Hawking temperature 

fc B W = ^. (10.3) 



- 25 - 



As we move along (that is, as w* increases, possibly even with w* — > +00), this 
relation continues to hold, with the Hawking temperature controlled by k(u*), as 
long as the adiabatic condition continues to hold. 

Working with this formalism we have never had to use any statements about 
any horizon forming anywhere — neither event horizon nor apparent horizon nor 
trapping horizon nor dynamical horizon — so this is consistent with all relevant 
previously known results § & g f§ § |, ft § |, [TJ, Q 0, © §§ ^ 0, [5T 



|32| , p3| , p4j |35| , pq , p7| , |40|| . For example, within this formalism one can check that the 
standard picture of Hawking evaporation is fully justified, even without making any 
particular commitment (pro or con) as to event-horizon formation. As the black hole 
loses mass its "surface gravity" k rises, and so does its Hawking temperature. And 
this picture should continue to hold all the way to the Planck scale. (At which stage 
it is the adiabatic condition that breaks down, along with the related breakdown 
of validity of the exponential approximation — "merely" indicating the need for a 
more thorough fully dynamic analysis.) The current formalism in no way "solves" 
the trans-Planckian problem and if anything brings it more sharply into focus. Over 
the years there have been repeated suggestions that Hawking emission should in some 
sense be "localized" (or more precisely "delocalized" ) to within an e-folding distance 
or two of the horizon (or in our case the might-be horizon). For the current state 



of affairs see [|53| p4| . Note that we have carefully described the Hawking flux as 
Planckian rather than thermal. To claim thermality one has to explicitly assume the 
formation of an event horizon (behind which one can hide correlations). 

Finally, it is interesting to remark that within this formalism other non-standard 
evaporation scenarios can also be envisaged and analyzed. For example, the evap- 
oration of hypothetical compact horizonless objects follows the described pattern, 
with a k(u) differing at each stage from ^Hawking HI- In addition, there is a 



chance that in these scenarios the evaporation process can proceed all the way to 
its end-point in a completely adiabatic manner |28|, 29|. Work on these important 
points is ongoing. 
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A. Surface gravity: peeling versus inaffinity 

In general, the peeling notion of surface gravity K pe eiing is defined by the "peeling off" 
properties of null geodesies near the horizon, while the more usual textbook notion of 
surface gravity /tmafimty is defined by the "inaffinity" of the naturally normalised null 
geodesic at the horizon. For Killing horizons these concepts coincide. For dynamical/ 
trapping/ apparent/ evolving/ putative horizons they will in general differ. In this 
appendix we shall motivate the relevant definitions — we will avoid the use of null 
coordinates, and adopt the more usual Schwarzschild curvature coordinates. 
For the static case, write the metric in the form [^5] 

ds 2 = -e- 2 * w [l - 2m(r)/r)]dt 2 + — — - + r 2 {d6 2 + sin 2 fl dcf 2 }. (A.l) 

1 — 2m(r)/r 

The Killing horizon is defined by the location where 2m(r)/r = 1, that is 

2m(r H ) = r H . (A.2) 



Then it is an old result [g6| that 

e -*"(l-2m' H ) 



2r H 

A.l The peeling notion of surface gravity 



(A.3) 



For the evolving case, write the metric in the form [^, [56|, |57 



ds 2 = -e- 2 ^ [1 - 2m(r, t)/r)]dt 2 + - — — + r 2 {d0 2 + sin 2 9 d<j) 2 }, (A.4) 

1 — 2m{r, t)/r 

and define the "evolving horizon" by the location where 2m(r,t)/r = 1, that is 

2m(r H (t),t) = r H (t). (A.5) 
To calculate K pe eimg note that a radial null geodesic satisfies 

(^y = e- 2 * w) [l-2m(r,t)/r)] 2 , (A.6) 

that is 



= ±e~^[l - 2m(r,t)/r)}. 
Now if the geodesic is near r#(t), that is r r^(t), we can Taylor expand 
dr e -^(r H (t),t)^ _ 2m'(r H (t),t)} 



(A.7) 



iJf - rn(f] [r(t)-r H (t)} + 0{[r(t)-r H (t)} 2 ), (A.8) 
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where now a prime denotes derivative with respect to r. That is, defining 

rCpeelingV'7 — 2rjf(t) ' V J / 

we have 

— = ±2 Kpeeling (t) [r(t) - r H (t)] + O ([r(t) - r H (t)} 2 ) . (A. 10) 

Then for two null geodesies ri(t) and r 2 (t) close to and on the same side of the 
evolving horizon 

d ^ Tl ~ T2) ~ ±2 Kpeeling (t) [n(t) - r 2 {t)}, (A. 11) 

which it is better to write as 



d\r 1 - r 2 \ 



2Kpeeling(t) \n (t) ~ T 2 (t) \ , (A.12) 



dt 

because this automatically keeps track of all the signs. Therefore 



\ri{t) - r 2 (t)\ w |ri(t ) - ^2(^0) I exp 



2 J ^peeling {t) dt 



(A.13) 



So 

^peeling (t) is a natural time-dependent generalization of the static version of Kpeeiing 
which correctly encodes the exponential "peeling off" behaviour of null geodesies. It 
is this quantity that we have seen is ultimately connected to the temperature of the 
Hawking flux. 

A. 2 The inafRnity notion of surface gravity 

In contrast, to calculate Kmaffinity consider the outward-pointing radial null vector 
field 

T= (l,e~* {rit) (l -2m(r,t)/r),0,0) . (A. 14) 

In a static spacetime, this null vector field is simply related to the Killing vector K a 
as 

r = K a + e a b K\ (A.15) 

where e a b is the 2-dimensional Levi-Civita tensor on the {r, t} plane. Hence the 
inaffmity «buik ( r ,t) defined by 

tVJ b = 2 «buik(r, t) £ b , (A.16) 

which always exists, everywhere throughout the spacetime, now naturally defines a 
notion of generalized surface gravity even for a time-dependent geometry. 



A brief calculation |)5|, |56|, |57[ shows that 



. . e \m(r, t)/r — m'fr, t)] & x . , . r „ „ , . . . 1 ■ , . , . 
Kbuik(r,t) = 1 1 ;/ r - e -*$'( r ,t)[l-2m(r,t)/r]--$(r,t). (A.17) 
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Note that at the evolving horizon 

, / \ \ e- q,(r »^[l-2m'(r H (t),t)} , . . 

K b ute(r H (t),t) = l — 1 HK h n - -$(r H (t),t) 

= Seeling(i) - -$(r H (t),t). (A.18) 

But this is not yet ftmaffinity (<0 • The quantity Kmaffimty(^) is found by somehow locating 
the true event horizon rg(t) 7^ (assuming an event horizon forms) and defining 

«inaffinity(0 = K bu i k (r B (t) , t) . (A.19) 

While we do not know exactly where the event horizon is, we can hope that when 
asymptotically approaching a quasi-static situation we have 

r E (t)^r H (t), (A.20) 

in which case we can expand in a Taylor series 

KWfinityW ~ «buik(r# (f) , t) + n' hulk (r H (t) , t) [r E (t) - r H (t)]. (A.21) 

That is 

/Wfinity(f) ~ «peeling(*) - ^ (r H (*) , t) + < u lk(^ (*),*) [n5 (*) - Hf(*)]- (A.22) 

This clearly demonstrates that the two concepts are in general distinct, though they 
will automatically agree in static (or stationary) spacetimes. Calculating the quantity 
K^ ulk (r//(t), t) is do- able but unedifying. 
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